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Abstract 

The cosmological gamma-ray burst (GRB) 
model for the production of uhra- high en- 
ergy cosmic rays is described, and the current 
observational evidence which support it dis- 
cussed. Several predictions of the model are 
presented, which would allow it to be tested by 
future high energy cosmic ray and gamma-ray 
experiments. If the predicted signatures of the 
GRB model are observed, they will not only 
corroborate the model, but will also provide 
information about the source population, and 
will allow to investigate the unknown structure 
of the inter-galactic magnetic field. 



1 Introduction 

Recent gamma ray and cosmic ray observa- 
tions give increasing evidence that the sources 
of gamma ray bursts (GRBs) and of cosmic 
rays with energy E > lO^^eV are cosmologi- 
cal (see [|| for GRBs; [2-5] for cosmic rays). 
The sources of both phenomena, however, re- 
main unknown. In particular, most of the 
cosmic ray sources discussed so far have dif- 
ficulties in accelerating cosmic rays up to the 
highest observed energies Furthermore, 
the arrival directions of the few cosmic rays 
detected above 10^°eV are inconsistent with 



the position of any astrophysical object that 
is likely to produce high energy particles 0, 
since the distance traveled by such particles 
must be smaller than lOOMpc ||]. 

Although the source of GRBs is un- 
known, their observational characteristics im- 
pose strong constraints on the physical condi- 
tions in the 7-ray emitting region which 
imply that protons may be accelerated by 
Fermi's mechanism in this region to energies 
1020 - lO^ieV 0, HI. The observed energy 
spectrum of cosmic rays above lO^^eV (UHE- 
CRs) is consistent with a cosmological distri- 
bution of sources of protons with a power-law 
generation spectrum typical of Fermi acceler- 
ation |]. Further more, the average rate (over 
volume and time) at which energy is emitted 
as 7-rays by GRBs and in UHECRs in the cos- 
mological scenario is, remarkably, comparable 
]To| , 1^. These facts suggest that GRBs and 
UHECRs have a common origin. 

We describe the GRB model for UHECRs in 
Sec. UJ], and discuss the current observational 



evidence which support it in Sec. 2.2. Several 



predictions of the model are presented In Sec. 
^. Sec. ^ contains a discussion. 



2 The GRB model 



1 



2.1 Fermi acceleration in dissipa- 
tive wind models of GRB's 

Whatever the ultimate source of GRBs is, ob- 
servations strongly suggest the following sim- 
ilar scenario for the production of the bursts 
The rapid rise time, ~ ms, observed in 
some bursts implies that the sources are com- 
pact, with a linear scale rp ~ lO^cm. The high 
luminosity required for cosmological bursts, 
L ~ lO^^ergs^^, then results in an initially op- 
tically thick (to pair creation) plasma of pho- 
tons, electrons and positrons, which expands 
and accelerates to relativistic velocities. In 
fact, the hardness of the observed photon spec- 
tra, which extends to ~ lOOMeV, implies that 
the 7-ray emitting region must be moving rel- 
ativistically, with a Lorentz factor 7 of order 
a few hundreds. 

If the observed radiation is due to photons 
escaping the expanding "wind" as it becomes 
optically thin, two problems arise. First, the 
photon spectrum is quasi-thermal, in contrast 
with observations. Second, the plasma is ex- 
pected to be "loaded" with baryons which may 
be injected with the radiation or present in the 
atmosphere surrounding the source. A small 
baryonic load, > lO^^M©, increases the opti- 
cal depth (due to Thomson scattering) so that 
most of the radiation energy is converted to 
kinetic energy of the relativistically expand- 
ing baryons before the plasma becomes opti- 
cally thin. To overcome both problems, it was 
suggested that the observed burst is produced 
once the kinetic energy of the ultra-relativistic 
ejecta is dissipated, due to collision of the rela- 
tivistic baryons with the inter-stellar medium 
or due to internal collisions within the ejecta, 
at large radius r = Vd > lO^^cm beyond the 
Thomson photosphere, and then radiated as 
7-rays @. 

Since 7^1, substantial dissipation of ki- 
netic energy at r = implies that the ran- 
dom motions in the wind rest frame are (at 



least mildly) relativistic. The relativistic ran- 
dom motions are likely to give rise to a tur- 
bulent build up of magnetic fields, and there- 
fore to Fermi acceleration of charged particles. 
We derive below the constraints that should be 
satisfied by the wind parameters in order to al- 
low acceleration of protons at the dissipation 
region to ~ 10^°eV. 

The most restrictive requirement, which 
rules out the possibility of accelerating par- 
ticles to energies ~ lO'^^eV in most astrophys- 
ical objects, is that the particle Larmor radius 
Rl should be smaller than the system size . 
In our scenario we must apply a more strin- 
gent requirement. Due to the wind expansion 
the internal energy is decreasing and therefore 
available for proton acceleration (as well as for 
7-ray production) only over a comoving time 
td ~ fd/^c. The typical Fermi acceleration 
time is ta — Rl/c leading to the require- 
ment Rl < Td/j. This condition sets a lower 
limit for the required comoving magnetic field 
strength 
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where E = WE2oeV, 7 = 3OO7300, L = 
lO^^Lsierg s^^ is the wind luminosity, and 
Be.p. is the equipartition field, i.e. a field with 
comoving energy density similar to that asso- 
ciated with the random energy of the baryons. 

The accelerated proton energy is also lim- 
ited by energy loss due to synchrotron radia- 
tion. The condition that the synchrotron loss 
time should be smaller than the acceleration 
time is 

B<3x 10^|oo^2o'G. (2) 

Since the equipartition field is inversely pro- 
portional to the radius r, this condition may 
be satisfied simultaneously with (|l|) provided 
that the dissipation radius is large enough, i.e. 
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The high energy protons lose energy also in 
interaction with the wind photons (mainly 
through pion production). It can be shown, 
however, that this energy loss is less impor- 
tant than the synchrotron energy loss pO[ |. 

The conditions imply, that a dissi- 

pative ultra-relativistic wind, with luminosity 
and bulk Lorentz factor implied by GRB ob- 
servations, satisfies the constraints necessary 
to allow the acceleration of protons to ener- 
gies ~ 10^'^eV by second order Fermi accel- 
eration, provided that turbulent build up of 
magnetic fields during dissipation gives rise 
to fields which are close to equipartition. It 
should be noted that equipartition field is also 
required in most of the dissipative relativis- 
tic wind models for GRBs. Gamma-rays are 
produced in these models by synchrotron and 
synchro-self Compton radiation of relativistic 
electrons produced by the dissipation shocks. 
Equipartition field is required in this case to 
allow efficient radiation of the electrons 

Although the details of the mechanism of 
turbulent build up of an equipartition field are 
not fully understood, it seems to operate in 
a variety of astrophysical systems. We note 
that a magnetic field of this strength may ex- 
ist in the plasma prior to the onset of internal 
collisions if a substantial part of the wind lu- 
minosity is initially, i.e. at r ~ ro, provided 
by magnetic field energy density, and if the 
field becomes predominantly transverse. The 
pre-existing field may suppress small scale tur- 
bulent motions. In this case shocks coherent 
over a scale R ^ rd/^ would exist and protons 
would be accelerated by first rather than sec- 
ond order Fermi mechanism. The constraints 
(0^) are valid in this case too, therefore leav- 
ing the above conclusions unchanged. 

2.2 UHECR spectrum and flux 

In the GRB model for UHECR production de- 
scribed above, the high energy cosmic rays are 



protons accelerated by Fermi's mechanism in 
sources that are distributed throughout the 
universe. In Fig. ^ we compare the UHECR 
spectrum, reported by the Fly's Eye and the 
AGASA experiments with that ex- 

pected from a homogeneous cosmological dis- 
tribution of sources, each generating a power 
law differential spectrum of high energy pro- 
tons dN/dE (X E'"^-^, as typically expected 
from Fermi acceleration (e.g. j^). (For 
this calculation we have used a flat universe 
with zero cosmological constant and Hq = 
75km s"""^; The spectrum is insensitive to the 
cosmological parameters and to source evolu- 
tion, since most of the cosmic rays arrive from 
distances < 500Mpc). The AGASA fiux at 
3 - 10 X lO^^eV is 1.7 times higher than that 
reported by the Fly's Eye, corresponding to a 
systematic ^ 20% larger estimate of event en- 
ergies in the AGASA experiment compared to 
the Fly's Eye experiment (see also [p| ^). We 
have therefore multiplied in Fig. npthe Fly's 
Eye energy by 1.1 and the AGASA energy by 
0.9. Bird et al. g find that the Fly's Eye fiux 
in the energy range 4x 10^^— 4 x lO^^eV can be 
fitted by a sum of two power laws: A steeper 
Galactic component with J cx E~^-^ domi- 
nating at lower energy, and a shallower extra- 
Galactic component with J cx E^^'^ dominat- 
ing at higher energy. The Bird et al. fit to 
the extra-Galactic component is also shown in 
Fig. I 

The data are consistent with the cosmolog- 
ical model for E > 2 x lO^^eV. Furthermore, 
the flux predicted by the model at lower en- 
ergy is consistent with the Bird et al. fit to the 
extra-Galactic component. (The flux deduced 
from the highest energy event in the Fly's Eye 
data is significantly higher than that predicted 
from the cosmological model. However, the 
statistical significance of the apparent discrep- 
ancy is not high: For the Fly's Eye exposure, 
the model predicts an average of ~ 1.3 events 
above 10^°eV, and the probability that the 
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Figure 1: The UHECR flux expected in a cos- 
mological model, compared to the Fly's Eye 
and AGASA data. Integers indicate the num- 
ber of events observed, la energy error bars 
are shown for the highest energy events. The 
dashed line denotes the fit by Bird et al. 
for the extra-galactic flux. 

first event observed at this energy range is 
above 2 x lO^OeV is ~ 15%.) The deficit in the 
number of events detected above 5 x lO^^eV, 
compared to a power-law extrapolation of the 
flux at lower energy, is consistent with that ex- 
pected due to a cosmological "black-body cut- 
off" . However, with current data the "cutoff" 
is detected with only 2a significance. 

The present rate at which energy should 
be produced as 10^^-10^°eV protons by the 
cosmological cosmic ray sources in order 
to produce the observed fiux is 4 ± 2 x 
10'*"'erg Mpc^'^yr"-'^. This rate is comparable 
to that produced in 7-rays by cosmological 
GRBs: The rate of cosmological GRB events is 
i/^ ~ 3 X 10"^Mpc~^yr"i Q, each producing 
« 3 X lO^^erg, corresponding to a 7-ray energy 
production rate of ^ lO'^^erg Mpc~^yr~^. 



The above analysis implies that the GRB 
model of UHECR production would produce 
UHECR fiux consistent with that observed, 
provided the efficiency with which the wind ki- 
netic energy is converted to 7-ray and UHECR 
energy is similar. There is, however, one addi- 
tional point which requires consideration. The 
energy of the most energetic cosmic ray de- 
tected by the Fly's Eye experiment is in ex- 
cess of 2 X 10^°eV, and that of the most en- 
ergetic AGASA event is above lO^^eV. On 
a cosmological scale, the distance traveled by 
such energetic particles is small: < lOOMpc 
(50Mpc) for the AGASA (Fly's Eye) event 
(e.g., 1^). Thus, the detection of these events 
over a ~ Syr period can be reconciled with the 
rate of nearby GRBs, ~ 1 per 50 yr in the field 
of view of the CR experiments out to lOOMpc 
in a standard cosmological scenario [Q, only 
if there is a large dispersion, > 50yr, in the 
arrival time of protons produced in a single 
burst (This implies that if a direct correlation 
between high energy CR events and GRBs, 
as recently suggested in ||l5|, is observed on a 
~ lOyr time scale, it would be strong evidence 
against a cosmological GRB hypothesis) . 

The required dispersion is likely to occur 
due to the combined effects of deflection by 
random magnetic fields and energy dispersion 
of the particles. Consider a proton of en- 
ergy E propagating through a magnetic field 
of strength B and correlation length A. As 
it travels a distance A, the proton is typi- 
cally deflected by an angle a X/Rl, where 
= E/eB is the Larmor radius. The typi- 
cal deflection angle for propagation over a dis- 
tance d is 9s ^ {d/\y^^\/RL. This deflection 
results in a time delay, compared to propaga- 
tion along a straight line, of order t{E, d) w 
Ogd/c Ri {eBd/E)'^X/c. The random energy 
loss UHECRs suffer as they propagate, owing 
to the production of pions, implies that at any 
distance from the observer there is some finite 
spread in the energies of UHECRs that are ob- 
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served with a given fixed energy. For protons 
with energies > lO^^eV the fractional RMS en- 
ergy spread is of order unity over propagation 
distances in the range 10— lOOMpc (e.g. Q). 
Since the time delay is sensitive to the particle 
energy, this implies that the spread in arrival 
time of UHECRs with given observed energy is 
comparable to the average time delay at that 
energy T{E,d). The field required to produce 
r > lOOyr is 

bVx > 10-"£;2orfrooG Mpc^/^ (4) 

where d = lOOdiooMpc. The required field 
is consistent with the current upper limit for 
the inter-galactic magnetic field, BX^^^ < 
lO-^G Mpc^/2 [[|. A time broadening over 
T 3> lOOyr is therefore possible. 

It should be noted, that a GRB producing 
3 X lO^^erg in UHECRs at 50Mpc distance, 
would produce a total fluence at Earth of ^ 2 
cosmic rays above lO^^eV per km^. In the 
presence of a magnetic field induced time de- 
lay, the typical distance dm{E) to the bright- 
est source observed over an energy range IS.E 
around E, with ^EjE ~ 1, is the radius of a 
sphere within which the average time between 
bursts is equal to the characteristic time de- 
lay T{E,d^[E)\; i.e. 4^d^z.^T(£;, d„)/3 = 1. 
Thus, the brightest source distance is 

drr,{E) ~ 30 { ^\ Ell' Mpc, 

(5) 

and its flux 

is / ~ 0.1^'/'(BAi/2/10-"G Mpci/2)-2/5 
per Wyr 11. Here E = lOi^SigcV. 



3 Predictions 



3.1 The Number and Spectra of 
Bright Sources 

The initial proton energy, necessary to have an 
observed energy E, increases with source dis- 
tance due to propagation energy losses. The 
rapid increase of the initial energy after it ex- 
ceeds, due to electron-positron production, the 
threshold for pion production effectively re- 
sults in a cutoff distance, dc{E), beyond which 
sources do not contribute to the flux above E. 
Since dc{E) is a decreasing function of i?, for a 
given number density of sources there is a criti- 
cal energy Ec, above which only one source (on 
average) contributes to the flux. For burst- 
ing sources, Ec depends on the product of the 
burst rate v and the time delay. In the GRB 
model, the burst rate is given by the GRB rate 
V = v^, which is determined from the GRB 
flux distribution. The time delay depends on 
the unknown properties of the intergalactic 
magnetic fleld, r cx B^X{d/E)^. However, 
the rapid decrease of dc{E) with energy near 
10^°eV implies that Ec is not very sensitive to 
the unknown value of B^X. For the range al- 
lowed for the GRB model, lO^^G Mpc^/^ < 
< 10-^G Mpc^/^ (the lower limit deter- 
mined by (U), and the upper limit by Faraday 
rotation observations ||l6|), the allowed range 
of Ec is lO^OeV < Ec < 3 x lO^°eY. 

Fig. H presents the flux obtained in one real- 
ization of a Monte-Carlo simulation described 
in ref. of the total number of UHECRs 
received from GRBs at some fixed time. For 
each realization the positions (distances from 
Earth) and times at which cosmological GRBs 
occurred were randomly drawn, assuming an 
average rate i^^ = 2.3x lO^^Mpc^'^yr^^, an in- 
trinsic generation spectrum np(E) cx E~^dE, 
and £'c = 1.4 X lO'^^eV. Most of the realiza- 
tions gave an overall spectrum similar to that 
obtained in the realization of Fig. |^ when the 
brightest source of this realization (dominat- 
ing at lO^^eV) is not included. At i? < E^ 
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Figure 2: Results of a Monte-Carlo realiza- 
tion of the bursting sources model: Thick solid 
line- overall spectrum in the realization; Thin 
solid line- average spectrum, this curve also 
gives dc {E) ; Dotted hues- spectra of brightest 
sources at different energies. 

the number of sources contributing to the flux 
is very large, and the overall UHECR flux re- 
ceived at any given time is near the average 
(the average flux is that obtained when the 
UHECR emissivity is spatially uniform and 
time independent). At E > Ec, the flux will 
generally be much lower than the average, be- 
cause there will be no burst within a distance 
dc{E) having taken place sufficiently recently. 
There is, however, a significant probability to 
observe one source with a flux higher than the 
average. A source similar to the brightest one 
in Fig. ^ appears ~ 5% of the time. 

At any fixed time a given burst is observed 
in UHECRs only over a narrow range of en- 
ergy, because if a burst is currently observed 
at some energy E then UHECRs of much 
lower energy from this burst have not yet ar- 
rived, while higher energy UHECRs reached 



us mostly in the past. As mentioned above, 
for energies above the pion production thresh- 
old, E ^ 5 X f O^^eV, the dispersion in arrival 
times of UHECRs with fixed observed energy 
is comparable to the average delay at that en- 
ergy. This implies that the spectral width AE 
of the source at a given time is of order the 
average observed energy, AE ^ E. Thus, 
bursting UHECR sources should have nar- 
rowly peaked energy spectra, and the brightest 
sources should be different at different ener- 
gies. For steady state sources, on the other 
hand, the brightest source at high energies 
should also be the brightest one at low ener- 
gies, its fractional contribution to the overall 
flux decreasing to low energy only as dc{E)~^ . 

3.2 Spectra of Sources at E' < 4 x 

lO^^eV 

The detection of UHECRs above W°e\ im- 
ply that the brightest sources must lie at dis- 
tances smaller than lOOMpc. UHECRs with 
E < Ax f O^^eV from such bright sources will 
suffer energy loss only by pair production, be- 
cause at < 5 X f 0^^ eV the mean-free-path 
for pion production interaction (in which the 
fractional energy loss is ~ f 0%) is larger than 
fGpc. Furthermore, the energy loss due to 
pair production over fOOMpc propagation is 
only ~ 5%. 

In the case where the typical displacement 
of the UHECRs due to deflections by inter- 
galactic magnetic fields is much smaller than 
the correlation length, A ^ d6g{d,E) ~ 
d{d/X)^/^X/RL, all the UHECRs that arrive 
at the observer are essentially deflected by 
the same magnetic field structures, and the 
absence of random energy loss during propa- 
gation implies that all rays with a fixed ob- 
served energy would reach the observer with 
exactly the same direction and time delay. At 
a fixed time, therefore, the source would ap- 
pear mono-energetic and point-like. In reality, 
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energy loss due to pair production results in 
a finite but small spectral and angular width, 
AE/E^de/9s < 1% |9). 

In the case where the typical displacement 
of the UHECRs is much larger than the cor- 
relation length, A ^ d9s{d,E), the deflection 
of different UHECRs arriving at the observer 
are essentially independent. Even in the ab- 
sence of any energy loss there are many paths 
from the source to the observer for UHECRs 
of fixed energy E that are emitted from the 
source at an angle 9 < 9s relative to the 
source-observer line of sight. Along each of 
the paths, UHECRs are deflected by inde- 
pendent magnetic field structures. Thus, the 
source angular size would be of order 9^ and 
the spread in arrival times would be compa- 
rable to the characteristic delay r, leading to 
AE / E 1 even when there are no random en- 
ergy losses. The observed spectral shape of a 
nearby (d < lOOMpc) bursting source of UHE- 
CRs at £^ < 4x lO^^eV was derived for the case 
A <C d9s (d, E) in , and is given by 



dN 

Ie 



n=l 



(-1) 



"+1 exp 



2n^n^E^ 



El{t,d) 



(6) 

where Eo{t,d) = de{2B'^\/ictY/'^ . For this 
spectrum, the ratio of the RMS UHECR en- 
ergy spread to the average energy is 30% 

Fig. ^ shows the line 9sd = \ the B — X 
plane, for a source at a distance d ~ 30Mpc 
observed at energy E ~ lO^^eV. Since the 
9sd = X line divides the allowed region in the 
plane at A ~ IMpc, measuring the spectral 
width of bright sources would allow to deter- 
mine if the field correlation length is much 
larger, much smaller, or comparable to IMpc. 

3.3 Correlation with Large Scale 
Structure (LSS) 

If the UHECR sources are indeed extra- 
Galactic, and if they are correlated with lu- 
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Figure 3: The line 9sd = A for a 
source at 30Mpc distance observed at en- 
ergy E ~ lO-'^^eV (dot-dash line), shown with 
the Faraday rotation upper limit i3A^/^ < 
10"^G Mpc^/^ (sohd hue), and with the lower 
hmit BA^/^ > 10-"G Mpc^/^ required in the 
GRB model. 



minous matter, then the inhomogeneity of the 
large scale galaxy distribution on scales < 
lOOMpc should be imprinted on the UHECR 
arrival directions. The expected anisotropy 
signature and its dependence on the relative 
bias between the UHECR sources and the 
galaxy distribution and on the (unknown) in- 
trinsic UHECR source density have been ex- 
amined in |2^. The galaxy distribution was 
derived from the IRAS 1.2 Jy redshift survey, 
which gives an acceptable description of the 
LSS out to - lOOMpc (see Q for a detailed 
analysis). 

Figures ^ presents a map of the angular de- 
pendence of the mean (over different realiza- 
tions of source distribution) UHECR intensity, 
for E > 6x lO^^eV and a model where UHECR 
sources trace galaxies with no bias. The 
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Figure 4: AitofF projection (Galactic coordinates) of the fractional deviation (from the 
all sky average) of the mean UHECR intensity. The heavy curve denotes the zero con- 
tour. Solid (dashed) contours denote positive (negative) fractional fluctuations at intervals 
[—0.5, —0.25, 0, 0.25, 0.50, 1.0, 1.5]. The super-galactic plane is denoted by the heavy sohd 
curve roughly perpendicular to the Galactic plane. The dotted curve denotes the Fly's Eye 
coverage of declination > —10°. The shaded regions show the high and low density regions 
used in the X{E) statistic. 



map clearly reflects the inhomogeneity of the 
large-scale galaxy distribution- the overdense 
UHECR regions lie in the directions of the 
"Great Attractor" [composed of the Hydra- 
Centaurus (/ = 300-360°, b = 0-+45°) and 
Pavo-Indus (/ = 320-360°, b = -45-0°) super- 
clusters] and the Perseus-Pisces super-cluster 
(/ = 120-160°, b = -30-1-30°). 

In order to determine the exposure re- 
quired to discriminate between isotropic 
and LSS correlated UHECR source distri- 
bution, the distribution of a statistic simi- 
lar in spirit to was considered, X{E) — 
J2i [ni{E) - nij{E)Y /nij{E), where ni is the 



number of events detected in angular bin I 
and nij is the average number expected for 
isotropic distribution (For the calculation of X 
24° X 24° bins were used; see also Fig. ^. From 
the X{E) distributions, generated by Monte- 
Carlo simulations of the UHECR source dis- 
tributions, it was found that the exposure re- 
quired for a northern hemisphere detector to 
discriminate between isotropic UHECR source 
distribution and an unbiased distribution that 
traces the LSS is approximately 10 times the 
current Fly's Eye exposure (0.1 the expected 
Auger exposure). If the UHECR source dis- 
tribution is strongly biased, in a way similar 
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to that of radio galaxies, the required expo- 
sure is ~ 3 times smaller. Furthermore, with 
10 times the current Fly's Eye exposure, it 
would be possible to discriminate between bi- 
ased and non-biased source distribution. The 
anisotropy signal is not sensitive to the cur- 
rently unknown number density of sources. 

Stanev et al. have recently reported 

that the arrival directions of > 4 x lO^^eV 
UHECR events detected by the Haverah Park 
experiment show a concentration in the direc- 
tion of the super-galactic plane. In agreement 
with Stanev et at, it is found in | pO[ that 
the probability to obtain the Haverah Park 
results assuming an isotropic source distribu- 
tion is very low. However, the results of 
show that this probability is not significantly 
higher for models where the source distribu- 
tion traces the LSS; thus, the concentration of 
the Haverah Park events towards the super- 
galactic plane can not be explained by the 
known LSS. It is important to note that for 
the biased model the probability to obtain the 
Haverah Park results is smaller than for the 
unbiased one. This reflects the fact that the 
super-clusters, while concentrated towards the 
super-galactic plane, have offsets above and 
below the super-galactic plane which cause the 
inferred UHECR distribution to be less flat- 
tened than seen in the Haverah Park data. 

4 Conclusions 

The GRB model for UHECR production has 
several predictions, which would allow it to be 
tested with future experiments. In this model, 
the average number of sources contributing 
to the flux decreases with energy much more 
rapidly than in the case where the UHECR 
sources are steady 0. A critical energy ex- 
ists, lO^^eV <Ec<5x lO^OeV, above which a 
few sources produce most of the flux, and the 
observed spectra of these sources is narrow, 



AE/E ~ 1: the bright sources at high energy 
should be absent in UHECRs of much lower 
energy, since particles take longer to arrive 
the lower their energy. In contrast, a model 
of steady sources predicts that the brightest 
sources at high energies should also be the 
brightest ones at low energies. 

At the highest energies, where most of the 
cosmic rays should come only from a few 
sources, bursting sources should be identified 
from only a small number of events from their 
coincident directions. Many more cosmic rays 
need to be detected at lower energies, where 
many sources contribute to the flux, in or- 
der to identify sources. Recently, the AGASA 
experiment reported the presence of 3 pairs 
of UHECRs with angular separations (within 
each pair) < 2.5°, roughly consistent with the 
measurement error, among a total of 36 UHE- 
CRs with E > 4 X IQi^cV The two 
highest energy AGASA events were in these 
pairs. Given the total solid angle observed 
by the experiment, ^ 27rsr, the probability 
to have found 3 pairs by chance is ~ 3%; 
and, given that three pairs were found, the 
probability that the two highest energy events 
are among the three pairs by chance is 2.4%. 
Therefore, this observation favors the bursting 
source model, although more data are needed 
to confirm it. 

Above Ec, there is a significant probabil- 
ity to observe one source with a flux consid- 
erably higher than average. If such a source 
is present, its narrow spectrum may produce 
a "gap" in the overall spectrum, as demon- 
strated in Fig. ||. It had recently been argued 
Ip^ that the observation of such an energy 
gap would imply that the sources of > lO'^^eV 
UHECRs are different from the sources at 
lower energy, hinting that these are produced 
by the decay of a new type of massive parti- 
cle. We see here that this is not the case when 
bursting sources are allowed, owing to the time 
variability. If such an energy gap is present. 
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our model predicts that most of the UHECRs 
above the gap should normally come from one 
source. If our model is correct, then the Fly's 
Eye event above 2 x 10^°eV suggests that we 
live at one of the times when a bright source is 
present at high energies. However, given the 
present scarcity of UHECRs, no solid conclu- 
sions can be drawn. With the projected Auger 
experiment |^ , the number of detected UHE- 
CRs would increase by a factor ~ 100. If Ec 
is 2 X lO^^eV, then a few bright sources above 
10^°eV should be identified. 

For the GRB model, the expected num- 
ber of events to be detected by the 5000km^ 
Auger detectors from individual bright sources 
at E lO^^eV is of order 100 (cf. eq. §); 
see also |19 ). The spectral width of these 



sources depends on the correlation length A 
of the inter-galactic magnetic field: Very nar- 
row spectrum, AE/E < 1%, is expected for 
A > IMpc, and a wider spectrum, AE/E ^ 1, 
is expected for A < IMpc (see Fig. ^. With 
energy resolution of ~ 10%, the Auger detec- 
tors would easily allow to determine the spec- 
tral width of the sources, and therefore to put 
interesting constraints on the unknown struc- 
ture of the magnetic field. 

If the distribution of UHECR sources traces 
the large scale structure (LSS) of luminous 
matter, large exposure detectors should clearly 
reveal anisotropy in the arrival direction distri- 
bution of UHECRs above 4 x lO^^eV. With 10 
times the current Fly's Eye exposure (0.1 the 
expected Auger exposure), it would be possi- 
ble to determine whether the sources are dis- 
tributed isotropically, or correlate with known 
LSS. Furthermore, it would be possible to de- 
termine whether or not the source distribution 
is highly biased compared to /i?j4 5 galaxies (as 
radio galaxies are). Thus, the anisotropy sig- 
nal would provide constraints on the source 
population. 

Finally, we would like to note another possi- 
ble signature of the GRB model, which was not 



discussed above. The energy lost by the UHE- 
CRs as they propagate and interact with the 
microwave background is transformed by cas- 
cading into secondary GeV-TeV photons. A 
significant fraction of these photons can arrive 
with delays much smaller than the UHECR de- 
lay if much of inter-galactic space is occupied 
by large-scale magnetic "voids" , regions of size 
> 5Mpc and field weaker than 10~^^G. Such 
voids might be expected, for example, in mod- 
els where a weak primordial field is amplified 
in shocked, turbulent regions of the intergalac- 
tic medium during the formation of large-scale 
structure. For a field strength ^ 4 x lO^^^G 
in the high field regions, the value required to 
account for observed galactic fields if the field 
were frozen in the protogalactic plasma, the 
delay of protons produced by a burst at a dis- 
tance of lOOMpc is ~ lOOyr, and the fluence of 
secondary photons above lOGeV on hour-day 
time scales is /(> E) ~ IQ-'^E^^^cm-^ ||. 
This fluence is close to the detection threshold 
of current high-energy 7-ray experiments. 
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